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Abstract 

We model the algorithmic task of geometric elimination (e.g., 
quantifier elimination in the elementary field theories of real and com- 
plex numbers) by means of certain constraint database queries, called 
geometric queries. As a particular case of such a geometric elimi- 
nation task, we consider sample point queries. We show exponential 
lower complexity bounds for evaluating geometric queries in the gen- 
eral and in the particular case of sample point queries. Although this 
paper is of theoretical nature, its aim is to explore the possibilities 
and (complexity-) limits of computer implemented query evaluation 
algorithms for Constraint Databases, based on the principles of the 
most advanced geometric elimination procedures and their implemen- 
tations, like, e.g., the software package "Kronecker" (see [H]). This 
paper is based on J7^ and is only a draft specially prepared for CESSI 
2006, representing work in progress of the authors. It is not aimed for 
publication in the present form. 



1 Introduction and summary 

The framework of constraint databases was introduced in 1990 by Kanel- 
lakis, Kuper and Revesz [17] as a generalization of the relational database 
model. Here, a database consists of a finite number of generalized relations, 
rather than classical relations. When we consider constraint databases over 
the real numbers, a generalized relation is finitely represented by a Boolean 
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combination of polynomial equalities and inequalities over the reals. These 
so called constraint formulas finitely represent possibly infinite sets in some 
real space R". Therefore, the constraint database model provides an elegant 
and powerful model for applications that deal with infinite sets of points in 
some finite dimensional real space, and is hence well-suited for modelling, 
e.g., spatial databases. For example, the spatial relation consisting of the set 
of points on the northern hemisphere together with the points on the equator 
of the unit sphere in the three-dimensional space can be represented by 
the constraint formula + y'^ + = 1 /\ z > 

The constraint model has been extensively studied by now and various 
logic-based query languages have been considered [IS] • First-order logic over 
the reals, F0(+, x , =, <, 0, 1), augmented with relation names to address 
(generalized) relations in the input database, is the standard query language 
for constraint databases. As an example we may consider, for an input 
database, which contains just one ternary relation represented by S, the 
query expressed by the F0(+, x, =, <, 0, 1, S')-sentence 

(3r)(Va;)(V?/)(V2)(5(x, y, z) ^ + y"^ + z'^ < r^). 

This query expresses that the three-dimensional spatial relation 5" is bounded. 

The standard way to evaluate this query on a particular database, e.g., 
the hemisphere above, consists in replacing the subexpression S{x, y, z) in 
the query expression by the formula x"^ + y'^ + z^ = 1 /\ z > Q and next to 
eliminate from the resulting F0(+, x, =, <, 0, l)-formula the quantifiers that 
were introduced by the query expression. In our example, this adds up to 
eliminating the quantifiers from 

(3r)(Vx)(V2/)(Vz)((x2 + y"^ + = I ^ z >Q) ^ x^ + y'^ + z^ < r"^), 

which would result in the value true. 

In this paper, we do not consider the traditional application domains of 
constraint databases, such as spatial databases, but focus on a completely 
different domain, namely geometric elimination theory. Hereto, we extend 
the constraint database model in the sense that we allow databases also to 
contain functions, rather than only relations. In this paper, a typical (input) 
database schema will be of the form 

(i?i, Rr] Fi, Fg) 

with relation names Ri {i = l,...,r) and function names Fj {j = l,...,s). 
The relation names are interpreted, following the context, by algebraic or 
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semi-algebraic sets and the function names by polynomial, or exceptionally, 
rational functions defined over the complex or real numbers. 

The reason to include function symbols is two-fold. Firstly, functions 
appear naturally as byproducts of quantifier-elimination procedures and it 
is therefore suitable to consider them as in- or outputs of such algorithms. 
So, it is natural to include them if we want to model appropriately the new 
application domain of geometric elimination theory. Typical examples of such 
functions are the determinant and the resultants of systems of n homogeneous 
equations in n unknowns in the linear and non-linear case, respectively. 

The second reason is based on a complexity argument and explains why 
it does not suffice to represent a k-ary function just by a. {k+ l)-ary relation 
that stores the graph of the function. When we extend the constraint data- 
base model with functions and likewise extend the first-order query language 
including function symbols for the representation of input functions, we are 
sometimes able to write queries more economically with respect to the num- 
ber of quantifiers. This leads in turn to more efficient evaluation of these 
queries. As an example consider a input schema containing a unary function 
symbol F. The first-order query formula 

y = F{F{x)) 

of F0(-|-, X, =, <, 0, 1, F) defines, for each interpretation of the symbol F by 
a unary real- valued function /, all tuples (x, y) of real numbers satisfying 
y — P{x). On the other hand, if we model the function symbol F by means 

of the graph of /, i.e., using a database schema containing a binary relation 
symbol i?, then the relation y = P{x) becomes first-order expressible by the 
FO(-F, X, =, <, 0, 1, i?)-query 

{3z){R{x,z)^R{z,y)), 

which contains a quantifier. Observe that the evaluation of this query cannot 
be done directly, it requires the elimination of this quantifier. 

First-order logic over the real or complex numbers extended with relation 
and function names, used to address input relations and functions, allows to 
define output relations in the traditional way. On the other hand, the creation 
of output functions requires an extension of first-order logic by special terms 
in order to specify these functions. To illustrate this, we consider the following 
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example. Let be given the F0(+, x , =, <, 0, 1, Fn, Fnn, F)-formula 



n 



{3xi) ■ ■ ■ {3Xn){/\Fii{ui, ... 



Um)Xl H h Fin{Ui, Un)Xn = 



(t) 



i=l 



n 



A\J XiJ^O) ^ F{ui, ...,Um) =0 



i=l 



in which the the input schema is given by the function names Fij = 
l,...,n), representing polynomial input functions fij{ui, ...,Um) and where 
the output schema is given by the term F representing a polynomial function 
f{ui, ...,Um)- The formula (f) may be interpreted as a specification of the 
function f{ui, ...^Um) by the requirement that the condition F{ui, ...,Um) = 
in (f) reflects that the linear homogeneous equation system given by the 
matrix ^{ui, ...,Um) = {fij{ui, ■■■,Um))i<i,j<n has a non-zero solution. An 
example function that satisfies this specification is the determinant of the 
matrix $. Another example of a function that satisfies the specification (f) 
is the square of the determinant. In Section ??, we will discuss generalizations 
of this example to non-linear systems of homogeneous equations and it turns 
out that possible interpretations of the output functions are resultants. 

We remark that the variables ui,...,Um and xi,...,Xn play a different 
role in the expression (f). Therefore, we shall use the following terminol- 
ogy: we call Ml,..., Mm parameters and xi,...,Xn variables. The idea be- 
hind this distinction is the following. Suppose we are given a concrete data- 
base instance (/ii(mi, ...,Mm), ...,/nn(Mi, ■■■,Um)) of the the database schema 
(Fii(mi, ...,Mm), ...,F„„(mi, ...,Mm))|^ This iustauce gives rise to a new data- 
base containing a single (m + n)-ary relation described by the expression 



This relation can be viewed as an algebraic family, which is parameterized 
by Ml, Mm, of database instances consisting of a single n-ary relation in the 
variables 

^Or better?: Suppose we are given a concrete algebraic family of database instances 
Um), /rm('"i, w„i)), parameterized by ui,...,it„i, over the the database 
schema (-Fn, i^„„). This way we avoid using parameterized schemas, which are lit- 
tle bit strange. See also next page. 
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The expression (f) is exemplary for the way we shall specify output func- 
tions in this paper. Thus, output functions will typically take parameters as 

arguments. 

We also remark that specifications like (f) are no longer first-order for- 
mulas, but rather higher-order expressions. On the other hand, we may view 
these expressions as queries in our new constraint database model (extended 
with function symbols). Indeed, a (constraint database) query is usually 
defined as a partial computable function that transforms a given * con- 
straint database instance (over a previously fixed input schema) into a new 
constraint database instance (over an also previously fixed output schema). 
Observe that the example described by the expression (f) fits in this setup. 
Indeed, for a given instance of the input symbols Fij = l,...,n) any 
instantiation of the function symbol F, satisfying (f), can be considered as 
an output of the corresponding query. 

In the sequel, any algorithm that produces an output satisfying a speci- 
fication like (t) will be considered as a query evaluation algorithm. 

Once the constraint database model is extended with functions *, the 
constraint database formalism may be used to describe the algorithmic task 
of geometric elimination. Hereto, we introduce the notion of geometric query, 
which is based on the distinction between parameters and variables, consid- 
ered before. The idea is to apply queries, which are specified as above, to 
algebraically parameterized families of input databases. Roughly speaking, 
a geometric query is a transformation of algebraic families of constraint da- 
tabases that is parameterization independent. In this sense, the query which 
transforms an given input instance (/ii(tii, Um), fnn{ui, Um)) of the 
the previous database schema {Fii{ui, Um), Fnn{ui, Um)) into the de- 
terminant of the matrix $ = {fij{ui, Um))i<i,j<n, satisfies the specification 
(f ) and is geometric. 

A central contribution of this paper is the conclusion that in the largest 
possible sense, any known or thinkable geometric elimination procedure can 
be modeled as the evaluation of a suitable geometric query. This allows 
us now to use descriptive specifications in order to formulate basic tasks of 
geometric elimination theory. Note that the descriptive specifications are 
much more fiexible and problem adaptive than the traditional operative ones 
as, e.g., the task of the computation of a certain determinant, resultant or 
Chow form. In order to make the statement of this conclusion precise, we 
have to concretize the data structures which will be used to represent the 
databases which occur as inputs and outputs of geometric query evaluation 
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algorithms, i.e., of elimination procedures. 

In view of progress made in [] concerning the complexity of geometric 
elimination algorithms, we shall describe the relations and functions of a 
given database instance by the complexity model of essentially division-free 
boolean-arithmetic and arithmetic circuits, respectively (see Section ?? and [] 
for details on this model). By the way, we observe that the classical represen- 
tation of database instances by quantifier-free F0(+, x, =, <, 0, l)-formulas 
(or F0(+, X, =, 0, l)-formulas) is contained in this complexity model. 

With this conceptual tool at hand, we shall be able to certify the intrinsic 
non-polynomial character of geometric elimination (Theorem ??). However, 
this complexity result does not exclude that particular elimination tasks, as, 
e.g., the evaluation of certain resultants, would be solvable in polynomial 
time. Motivated by practical applications, we shall therefore consider the 
particular elimination task of evaluating sample point queries. 

The classical data model for constraint databases [18] does not support 
data exploration and local visualization. Indeed, a quantifier-free formula in 
disjunctive normal form, describing the output of a query, allows the answer- 
ing of, for instance, the membership question, but it is does not allow an 
easjQ exhibition of the output, by, e.g., the production of sample points, or, 
for low dimensions, a visualization of the output. To increase the tangibility 
of the output, we suggest considering a new type of query that produces 
sample points. Furthermore, it could be desirable to support an exploration 
of the neighborhood of a given sample point. Both aspects, namely find- 
ing a sample point and exploration of its neighborhood, encounter a simple 
expression in case of rationally parameterizable algebraic or semi-algebraic 
varieties. Therefore, we shall extend the concept of sample point query to 
queries that return rationally parameterized families of polynomial functions 
as output. Such queries will be called extended sample point queries. In 
Section ??, we shall prove that extended sample point queries, associated 
to first-order formulas containing a fixed number of quantifier alternations, 
cannot be evaluated in polynomial sequential time by so-called "branching- 
parsimonious algorithms". This lower bound result suggest that further re- 
search on the complexity of query evaluation in constraint database theory 
should be directed towards the identification of database and query classes 
that have a strongly improved complexity behavior. As a pattern for the de- 

^Remark: By definable choice [23j you can produce sample points in FO. The correctness 
of this satement is in the easy. 
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velopment of such a theory, we suggest a new type of ehmination algorithms 
which are based on the notion of system degree and use non-conventional 
data structures (see [2l[3llIlEl[10l[IIl[T3l[l5in6l|20lEI]). 

This paper introduces a number of new concepts for constraint database 
theory that sometimes require certain notions from algebraic complexity the- 
ory, algebraic geometry and commutative algebra. These notions can be 
found in standard textbooks, such as [6] (algebraic complexity theory), pQ 
(commutative algebra) and p2] (algebraic geometry). The reader only inter- 
ested in database issues may read this paper while skipping these technical 
details (and in particular the rather involved proof of Theorems ?? and ?? 
below) . 

The remainder of this paper is organized as follows. 

2 Preliminaries on the constraint database 
model 

We define the notions of constraint database schema and instance. In Section 
??, we will discuss representations of database instances. The definitions of 
constraint database schema and instance are generalizations of the traditional 
definitions of constraint databases [18] that also allow polynomial functions 
to be included in a database. 

We assume the existence of an infinite set of relation names and function 
names. 

Definition 1 An constraint database schema is a finite sequence 
(i?i, ...,Rr; Fi, ...,Fs) of relation names Ri {i = 1, ...,r) and function names 
Fj [j = 1, s), where r and s are integers. To each relation and function 
name Ri and Fj, natural numbers ar{Ri) and ar{Fj) are associated, called 
the arity of Ri and Fj, respectively (we remark that constants are modelled 
by function names of arity zero). □ 

Let S = (i?i, Rr] Fi, Fg) be a constraint database schema. Further 
on, we shall be interested in expressing queries in first-order logic over the real 
numbers extended with the relation and function names appearing in S. We 
shall write fO{+, x,=, <,0, 1,5) for F0(+, x,=, <,0, l,i?i, Fi, ...,Fs) 

and F0(+, x,=,0,l,5) for F0(+, x , =, 0, 1, i?^; Fi, F,) for these 
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first-order languages with relation symbols Ri,...,Rr and function symbols 
We denote by R and C the sets of the real and complex numbers. 

Definition 2 Let S = {Ri, Rr] Fi, Fg) be a constraint database 
schema. An constraint database instance over 5 is a finite sequence 
(Ai, Ar] fi, fs]snch that Ai is a semi-algebraic subset of R'^^(^>) or a 
constructible subseloof C"'''-'^*-' {i = 1, ...,r) and such that fj is a polynomial 
or rational function from R'"'(-^') to R or from C°^*^^'^ to C (j = 1, s). □ 

Below, we will refer to an constraint database schema and instance simply 
as database schema and database instance. 

Data objects such as the semi-algebraic and constructible sets and poly- 
nomial functions can be modeled in various ways. In the sequel, we shall use 
the term data model to refer to a conceptual model that is used to describe 
data objects. Quantifier-free first-order formulas over the reals represent an 
example of a data model which describes semi-algebraic sets. We use the 
term data structure to refer to the actual structure that implements the data 
models. For example, quantifier-free formulas may be given in disjunctive 
normal form and the polynomials appearing in them may be given in dense 
or sparse representation. 

It is important to remark that there are different data structures for 
constraint database instances and that for each particular structure there are 
different representations of particular instances. For example, relation and 
function instances can not only be finitely described by means of (sometimes 
unique) first-order formulas over the real numbers but also by, e.g., {never 
unique) essentially division-free arithmetic-boolean circuits and arithmetic 
circuits, respectively (see Section ?? and ??). For the purpose of this section, 
it is enough to assume some fixed data structure. When we speak about a 
representation of a database instance, it will be with respect to this fixed 
data structure. The reader should think that representations of database 
instances in this data structure are typically not unique. 

Next, we define the notion of query. 

Definition 3 Let an input schema S = {Ri, Rr] Fi, Fg) and an output 
schema S = {Ri, Rp] Fi, Fg) be given. 

■^Maybe Chevally's QE for algebraically closed fields should be mentioned. The con- 
straint database audience is probably not familiar with that. 
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A query over the input schema S and output schema 5 is a partial map- 
ping Q that maps database instances over S to database instances over S. 
This mapping Q can also be interpreted as a series of partial mappings 
{Qfij^, ■■■,Qjip', Qfj^, ■■■,QfJ, where Q^, {i = l,...,p) maps database instances 
over S to semi-algebraic or constructible subsets of R"''(^») or C"''^^'-' and 
where Qp. {j — 1, g) maps database instances over <S to a polynomial or 

rational function from R"''^^:') to R or from C"''^^^^ to C. We shall always sup- 
pose that queries are induced by partial mappings which map representations 
of input instances to representations of output instances. Hence, semantically 
equivalent representations of input instances (interpreted as semi-algebraic 
or constructible sets or polynomial or rational functions) are mapped to se- 
mantically equivalent representations of output instances. These mappings 
of representations of input instances or output instances always will be given 
implicitly by the context and we shall not specify them further. The reader 
may assume that these mappings are computable (in some suitable sense), 
but we will not rely on this fact in this paper. □ 

We use the terminology relational query and functional query in case of 
q = and p = 0, respectively. 

Below, we will simply refer to these mappings as queries whenever the 
input and output schemas are clear. 

We will be especially interested in relational and functional queries 
that are expressible in extensions of first-order logic over the real num- 
bers, F0(-|-, X , =, <, 0, 1), or complex number, F0(-|-, x , =, 0, 1). Indeed, 
given an input schema S = {Ri, Rr] Fi, Fs) and an output schema 
S = {Ri, Rp] Fi, Fq), we can consider F0(-|-, x, =, <, 0, 1, 5, 5) and 
F0(-|-, x,=,0,l,<S,<S) as logics to express relations between the input and 
output databases or as formalisms to specify relations and functions. In 
particular, any formula ip with k free variables in F0(+, x , =, <, 0, 1, 5) 
or F0(-|-, X , =, 0, 1, iS), when evaluated on a database instance over S = 
(i?i, Rr, Fi, Fg) defines a /c-ary output relation, when we interpret vari- 
ables to range over the real or complex numbers. 

The formula (p may be considered as a relational query that corresponds 
to an output schema S which consist of a single k-ary relation symbol. 
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3 Geometric elimination algorithms modeled 
as geometric queries 

In this section, we argue by means of a number of examples of elimination 
problems that geometric elimination algorithms can be modeled as geometric 
queries (to be defined further in this section) that satisfy some precise restric- 
tions. We also discuss the constraint formalism as a specification language. 

3.1 Elimination algorithms for non-parametric elimi- 
nation problems 

As a first example, let us consider the family (for varying polynomial function 
g) of elimination problems 

n 

{3xi)---{3Xn){/\x^i -Xi = Ag{Ui,...,Um,Xi,...,Xn) = 0), (t„) 
i=l 

in which ui, ...,Um are considered parameters. Elimination problems of the 
kind (ta) are typically solved by algorithms that produce as intermediate 

results suitable matrices Mj,-^, Mj,^ G Q^^^, which only depend on the 
system A^^-^^x'f — = (in case of the elimination problem (fa), we have 
A'^ = 2"). These matrices M^^^, Mx„ are defined in such a way that the 
parameters ui, ...,Um satisfy the first-order formula (fa) if and only if they 
satisfy the quantifier-free formula 

det^(Mi, ...,Um,Mr,^, ...,M^J = 0. (tb) 

By the way, let us observe that the formula (fb) represents a declarative 
specification of the matrices M^.^, M^.^. 

If the system A'^^-^xf — Xi — Ois changed by one that defines another zero- 
dimensional variety, in which polynomials of degree at most d appear, then 
N will also change and in fact be bounded by (T'. More generally, formulas 
of the form 

{3xi) ■ ■ ■ {3xn){ip{xi, Xn) A g{ui, Um, Xi, . . . , x^) = 0), (t) 

where (p{xi, ...,Xn) defines a zero-dimensional Q-definable subvariety of C", 
and g is a. polynomial function, form a well-known class of elimination prob- 
lems. Algorithms which solve these elimination problems produce typically 
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matrices M^.^, M^^ as output. These matrices depend in size and content 
only on the subformula ^p{xi, .... .t„) and they provide an easy way to obtain 
an ehmination formula, namely dei g{ui, ....,Um-, M.^^. M^^) = 0, for arbi- 
trary g. These algorithms can therefore be seen as a pre-processing of the 
formula (fi{xi, ...,Xn) into the matrices M^^, Mx„, aimed to facilitate the 
expression of the solution det g{ui, ...,Um, M^^^, M^^) = of the solution 
of the elimination problem (f). 

If we want to model this class of elimination problems and their elimina- 
tion polynomials using the language of constraint databases, there are several 
ways of doing this. Here, we start with one simple formulation and then we 
explain how the formalism of constraint databases can be flexibly adapted 
to deal with more general situations. 

For instance, we may produce the elimination problem (fa) by applying 
the query expressed by the following F0(-|-, x, =, 0, 1, Fi, Fn, (j)-formula 

n 

(3xi) ■ ■ ■ {3Xn)i/\Fi{xi, ...,Xn) = A G{ui, ...,Um,Xi, . . . , Xn) = 0) (fc) 
i=l 

over the input schema S = (Fi, G), to the database instance 
(/i, /„, g) given by fi{xi, x„) := xj - Xi {i = 1, n) and some polyno- 
mial g{ui,...,Um,Xi,...,Xn). 

Similarly, the more general elimination problem (f) may be produced by 
the query that is expressed by the F0(-|-, x, =, 0, 1, R, G)-formula 

{3Xi) ■ ■ ■ {3Xn){R{xi, Xn) A G{ui, Um,Xi,..., X^) = 0) 

over the input schema S' = {R,G), to the database instance {A,g), 
where A is supposed to be the zero- dimensional Q-definable subvariety of 
C" which is given by the F0(-|-, x, =, 0, l)-formula Xn) and where 

g{ui, Urn, xi, Xn) is some polynomial. In both cases, the output schema 
is iS = {A}, where A is a function symbol, with ar(A) = m. The database 
instance of S which represents the output of the query, is in both cases, the 
polynomial det m^, M^^, M^J. 

Another way to produce the elimination problem (f ) is the following: we 
suppose now that g{ui, Um, xi, Xn) is a fixed polynomial and we consider 
the query Qg expressed by the F0(-|-, x, =, 0, 1, i?)-formula 

(3xi) • • • {3Xn){R{Xi, Xn) A g{ui, Um,Xi,..., Xn) = 0) 
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over the input schema T = (R). 

We restrict our attention to input instances given by F0(+, x,=,0,l)- 
formulas f{xi, ...,x„) that define zero-dimensional varieties of fixed (strictly 
positive) cardinality N. Consequently, the output matrices Mx^, M^^ have 
also fixed dimension A^. In this way, we may view the output (M^jj, M^^) 
of the underlying elimination algorithm as an instance of a schema T = 
{Fi, Fi), with ar{Fj) = {j = 1, ...,£), where £ = n-N"^. Here, the function 
names Fj {j = !,...,£) are designed to describe the constant entries of the 
matrices Ma,^, M^^ (below, we will slightly abuse the notation by regarding 
{Mx^, Mx„) as an instance of the schema T = (-Fi, F^)). Observe that 
this defines a query Q with input schema T and output schema T which 
maps the given Q-definable input instance A C C" of cardinality N to 
QiA) := (Mx„...,MxJ. 

On the other hand, elimination problems like (f^), are obtained by apply- 
ing the queries Qg, expressed by quantifier-free F0(+, x , =, 0, 1, Fi, F^)- 
formulas 

dctg{ui,...,u„„Fi,...,Fe) =0 (U) 

over the input schema S = (Fi,...,^^) to the database instance containing 
the £ = n ■ entries of the matrices M^^, .... M^^^ as constants. 

Hence, the original input database A becomes pre-processed into a new 
databases consisting of the ?7.-tuple of matrices {M^-^, M^^) in such a way 
that, for arbitrary g^, the query Qg applied to the input database A and the 
query Qg applied to (M^.^, ...,Mr^) describe the same sets. For any input 
instance A, we have therefore 

Qg{A)^Qg{Q{A)). 

We observe that the pre-processing performed by the underlying elimination 
algorithm can be seen as the computation of a view on the original database 
instance A, that allows the replacement of the query Qg., which is defined us- 
ing quantifiers, by the query Q^, which can be expressed without quantifiers. 

Below, we shall pay particular attention to the following variant of the 
family of elimination problems (fa), namely 

n 

(3xi) • • • {3xn){/\ X- - = A y = g{ui, u^, Xi, . . . , (*„) 

i=l 

Here, Ui,...,u.m are considered as parameters, y as free variable and 
g{ui,... ) is a suitable polynomial function. Let M^^, e 
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q2"x2" ^YiQ matrices introduced before. Then all previous comments 
remain valid, mutatis mutandis, if one replaces the determinant in the 
quantifier-free formula (fh) by the characteristic polynomial (in the free vari- 
able y), of the matrix 

3.2 Elimination algorithms for parametric elimination 
problems 

In this section, we discuss three variations of the previous examples of elim- 
ination problem. 

3.2.1 First variation 

We are going to consider the following variation of the example (f) of Sec- 
tion I3.lt where the elimination problem takes the form 

n 

(3xi) • ■ ■ {3Xn)i/\Mui, ...,Um,Xi, ...,Xn) = A g{Ui, ...,Um,Xi, . . . , X„) = 0). 

1=1 

if) 

Here ui,...,Um are parameters, possibly subject to some first-order defin- 
able restriction (the parameters instances that satisfy this restriction are 
called admissible for the query (f)) and fi,---,fn,g are polynomials in 

Ul, . . . , Urm Xl, Xn- 

For fixed polynomial g{ui, ...,Um,Xi, ...,Xn) the elimination problem 
(f) can be produced by applying the query Q'g, expressed by the 
F0(+, x,=,0,l,Fi,...,F„)-formula 

n 

(3X1) ■ ■ ■ {3Xn){/\Fi{ui, ...,Um,Xi, ...,Xn) = A g{ui, ...,Um,Xi, . . . , X„) = 0) 

1=1 

(t'c) 

over the input schema S = (Fi,...,F„), where now ar{Fi) = m + n {i = 
1, r?,), to some suitable database instance (/i, ...,/n). 

For the sake of consistency and conciseness of exposition, we shall require 
that the input database should satisfy the following two flatness conditions: 

• the polynomials fi, fn form a regular sequence in 
Q [til , . . . , UjYi , Xl, Xji] , 
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• for V := {/i = 0,...,/„ = 0} C C™ x C", the morphism of affine 
varieties n : V ^ C", induced by the canonical projection C™" x C" — )■ 
C"', is finite. 

The ehmination problem (f) is typically solved by algorithms that pro- 
duce matrices with rational entries M^^^ui, ...,Um), M^^^ui, ...,Um) € 



We observe, that the flatness conditions satisfied by the input database 
(/i)---)/n) imply that the characteristic polynomial (and in particular the 
determinant) of the matrix g{ui, ...,Ufn,Xi, becomes a polynomial ex- 

pression in Ui, Urn- 



Therefore, if we wish to describe, in a similar way as in Section 13. 
such an elimination algorithm as a query mapping input database instances 
(/i) fn) to output database instances {M^^^ui, li^), M^^^ui, Um)), 
the constraint database model should be extended and allow an output 
schema that is not apriori fixed, e.g., using dynamic arrays of function names 
in this case. 

Since our actual constraint database model is limited to fixed schemas, 
we are obliged to require that the input database should satisfy a third flat- 
ness condition, namely that the typical fiber of the finite morphism of affine 
varieties n : V ^ C" should be of cardinality A^, where is a previously 
fixed, strictly positive integer. 

Then we may fix an output schema S = {Fi, Fi), as before, where 
i = n ■ and ar{Fj) = m {j = 1, i). An instance of S is now given by i 
rational functions belonging to Q{ui, ...,Um)- 

We call a parameter point {ai,...,am) € C"^ admissible if the (ra- 
tional) entries of M^^^^ui, ...,Um), Mx„{ui, ...,Um) are well-defined in 
am.) • Let us finally remark that the entries of the output matrices 
{Mj;^{ui, ...,Um), Mx„{ui, ...,Um)) are not arbitrary rational functions in 
Ui, ...jUn- Indeed, for every input instance (/i, fn) over S, and every two 
admissible parameter instances (cti, ...,«„) and {a[, ...,a^) for which fi{ai, 

am, Xi, Xn) and fi{a[, a'^, Xi, x„) {i = 1, n) are the same poly- 
nomials (in xi, Xn), also M^^^ai, am) and Mrc.{a[, a'm) will be equal. 
This remark will turn out to be crucial for the lower bound results that follow 
in the next section and motivate us to introduce the notion of "geometric 
query" at the end of this section. 



Q(mi, ...,Um) 



NxN 
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3.2.2 Second variation 

Let fi,...,fn be polynomials in the indeterminates ui, ...,Um,xi, ...,Xn and 
assume that are homogeneous of degrees di,...,dn with respect to 

The elimination problem 

n n 

) = A V X, 7^ 0) 

can be obtained by applying the query, expressed by the 
F0(+, X , =, 0, 1, Fi, F„)-formula 

n n 

{3x,)---{3xn){/\F,{ ) = A Y X, ^ 0) 

1=1 1=1 

over the input schema S = {Fi, ...,Fn} with ar{Fi) = m + n {i = 1, ...,n) 
to the database instance Then the output of an elimination 

algorithm could be the resultant of the polynomials with respect 

to Xi,...,Xn, which we denote by denoted by Res^^';;;'!^'^(/i, /„). Remark 
that parametric systems of n homogeneous linear equations in the unknowns 
xi, ...,Xn represent a particular case of this situation. In this case the resul- 
tant becomes the determinant. 

We observe now that for every input instance over S = 

(Fi,...,F„) and every two parameter instances {ai,...,am) and {a[, a'^) 
for which /^(ai, xi, x„) and fi{a[, a'^, Xi, Xn) (i = l,...,n) 
are the same polynomials, also Res^^';;;'f'^(/i(Q!i, a^, a;^), /„(«!, 

Qmj ^li •••) ^n)) and RgS^^' '.J," (/i(g;-|^, Cl^, Xi, X„), fniS^l-, ■•■1 O^mi 

...,Xn)) are the same polynomial. 

3.2.3 Third variation 

As a last variation on of the example of Section 13.1^ let us reconsider the 
family (for varying function g) of elimination problems 

n 

{3xi) ■ ■ ■ {3Xn){/\x'^^ -Xi = OA g{Ui, ...,Um,Xi, ...,Xn)= 0), (ta) 
i=l 

where the part A'^^^xj — = is considered as fixed. This example of a fam- 
ily of zero- dimensional elimination problems may seem trivial, but in fact. 
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no polynomial-time algorithm is known to solve it. It turns out that this ex- 
ample already illustrates the intrinsic difficulty of geometric elimination. For 
instance, elimination algorithms which use simultaneous Newton iteration to 
all zeroes of a zero-dimensional equation system lead to the consideration of 
this kind of problem, see, e.g., []. 

The complexity of the most efficient known elimination algorithms for 
problem (fa) depends exponentially on n, and linearly on the circuit com- 
plexity of g. 

The elimination problems (fa) can be obtained by applying the query 
expressed by the F0(+, x,=,0,l, G)-formula 

n 

{3xi) ■ ■ ■ {3Xn){/\x1 - Xi = A G{Ui, ...,Um,Xi, . . . , Xn) = O) (f^') 

i=l 

over the input schema S = (G), with ar{G) = m + n , to the data- 
base instance {g) consisting of the polynomial function g. The matrices 
Mj;^,...,Mj^ G Q^"^^", introduced in Section [STTl are well-defined in this 
case and the elimination algorithm under consideration produces as output 
a polynomial of the form det g{ui, ...,Um, M^^^, ...,Mx^). So, the elimination 
algorithm can be modeled by a query with input schema S = (G) and output 
schema S = (G), where ar{G) = m. 

Finally, we remark for any two admissible parameter instances 
(«!,..., dm) and {a[, a'^) for which g{ai, am, xi, Xn) and g{a[,..., 

the same polynomials (in xi, ...,Xn), also 
det^(Q;i, ...,am, M^^, M^J and det g{a[, a^, M^.^, M^J are the same 
values. 

3.3 Definition of geometric queries 

The obserations made earlier in Section [3] motivate the notion of geometric 
query, which we shall introduce below. First, we shall formalize a distinction 
between parameters and variables. 

3.3.1 Variables versus parameters 

From now on, we shall distinguish between two types of indeterminates which 
we refer to as variables on the one hand and as parameters on the other hand. 

The difference will be reflected in the notation: we use ui,U2, ■ ■ ■ to indi- 
cate parameters and Xi,X2, ■ ■ ■ to indicate variables. In this section, we shall 
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limit our attention to this purely syntactic distinction between variables and 
parameters and postpone the discussion of their (difficult) semantic meaning 
to Section m 

In particular, when dealing with geometric elimination, we shall always 
work with an input database schema S = {Ri, Rr', Fi, Fg) where all 
relation and function symbols are of arity m + n and which are such that 
the first m indeterminates refer to parameters and the last n to variables. 
Our notation reflects this circumstance by by . . . , Um', xi, . . . , x„) (^ = 

l,...,r) and . . . , w^; Xi, . . . , x„) {j = l,...,s). 

Let Ai be instances of Ri and /-,• be instances of Fj. We shall 
always assume that there are F0(+, x , =, <, 0, l)-formulas describing 
Ai{ui, . . . , Um', Xi, . . . , Xn) and polynomials (in the variables xi, x^) over 
the function field Q(ui, Um), describing fj{ui, . . . , Um] xi, . . . , x„) respec- 
tively. For an admissable instance (cti, of the parameters ui, ...,Um, 
we shall always require that /^(ai, cxm, xi, x„) is a well-defined polyno- 
mial in xi, Xn- 

Definition 4 Let S = {Ri, Rr] Fi, Fg) be a database schema and let 
(^1, Ar] fi, fs) be a database instance over S. 

We call two admissible instances (ai, . . . , am) and {a'l, . . . , a'^) of the pa- 
rameters Ui, ...,Um equivalent in the given database {Ai, Aj.; fi, fs) (or 
simply [Ai, A^; fi, fs)- equivalent) if the following conditions are satis- 
fied: 

• {{xi,...,Xn) e W I Ai{ai, . . . ,am,Xi, ...,Xn)} = {{Xi,...,Xn) G R'' | 

Ai{a'i, . . . , a'm, xi, Xn)} for all i = 1, r; and 

• fj{ai, . . . ,ani,xi, ...,Xn) and fj{a[,...,am,xi,...,Xn) are the same 
polynomials (in xi, Xn) for all j = 1, s. □ 

3.3.2 Definition of geometric queries 

Definition 5 Let S = (i?i, Rr] Fi, Fs) be an input database schema, in 
which all relations and functions are of arity m-|-n as discussed above, and let 
S = {Ri, Rp] Fi, Fg) be an output schema. For the sake of simphcity of 
exposition, we assume that the relations and functions appearing in instances 
over S and S depend on the same parameters ui, ■..,Um- 

A query Q = {Qj^_^, ...,Qj^^;Qp_^, ...,Qp^) over the input schema S and 

output schema S is called a geometric query if for any database instance 
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{Ai, Ar] /i, /s) and any two admissible (^i, Ar] fi, /s)-equivalent 
instances (ai, . . . , a^) and {a[, . . . , a^) for the parameters Ui, Um, the 
parameter instances (ai, . . . , am) and {a[, . . . , a'^) are also admissible and 
equivalent in the database 

Q{Ai,...,Ar;fi,...Js) = 

Q_Fl(^l) •••) ^r! /l) /s)) Q_Pg(^l) •••5 ^r! /l i •••) /«))• 

□ 

In the case of arbitrary first-order expressible queries, it may occur that 
the parameters appearing in the relations and functions in the output schema 
are different from those appearing in the relations and functions of the input 
schema. Nevertheless, this fact does not really restrict the applicability of our 
notion of geometric query, since the formula that defines the given first-order 
query may easily be rewritten (unnested) into an equivalent one in which the 
input and output schema depend on the same parameters. Since we will not 
relay on this more general situation, we shall not go into more details. 

It is clear that our notion of geometric query captures the reality of ex- 
isting elimination algorithms when the query-formula only uses symbols be- 
longing to the given database schema, as, e.g., (fc) in Section 13.11 If the 
query-formula contains both subexpressions without database symbols and 
subexpressions with database sysmbols, as, e.g., (f^) in Section 13.2.11 and 
(Yc) ill Section [3. 2. 3[ then the notion of geometric query still has an intuitive 
sense for elimination algorithms. In the general case, the query-formula must 
possibly be rewritten, before the formal definition of geometric query may 
be meaningfully applied. 

We end this section with some trivial examples of geometric queries. The 
transformation that maps the function U1U2X1 to the function xi is an exam- 
ple of a geometric query, as is the identity and the transformation that maps 
the functions UiXi,U2Xi to U1U2X1. 

The transformation that maps the function U1U2X1 to the function Ui is an 
example of a non-geometric query. Taking {ai,a2) := (1,0) and {a[,a2) = 
(0,1) gives two equal functions for what concerns U1U2X1, but not for for 
what concerns Ui. □ 
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3.4 Specification of queries 
3.4.1 First-order specification 

First-order logic over the real or complex numbers extended with relation 
and function names, used to address input relations and functions, allows to 
define output relations in the traditional way. 

As remarked in the introduction, the creation of output functions re- 
quires an extension of first-order logic by special terms in order to specify 
these functions. The F0(-|-, x , =, <, 0, 1, Fn, F)-formula (f) in the 
Introduction is a first example of such a specification. Both the determinant 
and the square of the determinant of any database instance of the schema 
(Fii, Fnn) satisfy this specification. 

A less trivial example is given by the specification formula 

n 

, Um, Xi, Xn) = 

1=1 

Fi{Ui, Urn, ZXi, ZXn) = 0) 

n 

= (t) 

i=l 

n 

a\/ Xi ^ 0)) ^ F{ui, Mm) = 
1=1 

over the input schema S = {Fi, ...,Fn), with ar{Fi) = m + n {i = 1, ...,n) 
and output schema S = (F) with ar{F) = m. 

This formula can be interpreted as a descriptive specification of a (typi- 
cally geometric) query which maps instances (/i, fn) over S = (Fi, F„) 
to instances (/) of the output schema S = (F). 

The reader should be aware that this formula, contrary to the formula 
(t) of the Introduction, does not apriori restrict the degree of the input 
polynomials with respect to the Xi, ...,Xn- 

However, fixing an input that satisfies the left-hand side 

of the equivalence in the specification (f), the polynomials fi satisfy for 
suitable non-negative integers di the condition fi{ui, ...,Um, zxi, zXn) = 
z'^'fi{ui, Urn, xi, Xn) {i = I, ...,n). As we have already observed in Sec- 
tion |32II1 for each occurrence of degrees {di, dn) the notion of resultant is 
well defined and denoted by Resj^'''''x^(/i, ■■■,fn)- Hence if we interpret the 
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function symbol F by Res^i'''''f" (/i, ■■■,fn), we obtain the output of a query 
which satisfies the given specification. 

The traditional informal specification of the resultant is operative, but as 
illustrated here, the formalism of constraint databases allows the declarative 
specification of objects such as resultants and determinants. 

3.4.2 Higher-order specification 

Now, we consider an example of a specification of an elimination task that 
cannot be formulated by means of a first-order formula. We use a specifica- 
tion in English, based on database schemas as before. 

As, we have already seen in Section 13.11 and I3.2.3[ for varying polynomial 
function g{ui, Um, Xi, ...,Xn), the family of elimination problems expressed 

by 

n 

(3xi) ■ • • i3Xn)i/\x'^i - Xi = Ay = g{ui, ...,Um,Xi, . . . (*a) 
i=l 

can be interpreted as the application of the query expressed by the 
F0(+, X, =, <, 0, 1, G)-formula 

n 

(3Xi) ■ ■ ■ {3Xn){/\xl -Xi = OAy = G{Ui, ...,Um,Xi, . . . (t'^') 

over the input schema S = {G) to the database instance (g) containing the 
polynomial function g. 

The canonical elimination polynomial of {*a) is 

p{y,ui,...,Um) ■■= Yl (y - g{ui,...,Um,ei,...,en)). (*e) 

(£i,...,en)G{0,l}" 

The elimination algorithms, described in the geometric elimination liter- 
ature, aim at obtaining a representation of the polynomial p, such that for 
any given values of the parameters Ui, Um, P can be evaluated as a polyno- 
mial function in y, e.g., by means of a division-free arithmetic circuit. More 
specifically, the aim is to find polynomials ui^ui, ...,Um), ■■■,0Je{ui, ...^Um) and 
a polynomial 

q{ti,...,ti,y) = gj{ti,...,te)y^, 

0<i<2" 
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belonging to Q[ti, ti] and Q[ti, t^, y] respectively, such that the identity 



p{y,Ui, ...,Um) = q{uJl{Ui, ...,Um), ...,Ui{Ui, ...,Um),y) 

2" 

= ■■■,Um), ■.■,iVe{Ui, ...,Um))y^ 

j=0 

is satisfied. 

The idea behind this elimination strategy is that of partial evaluation. 
This means that the evaluation of p with respect to y is postponed and 
that p{y,Ui, ...,Um) is written, for some fixed values of Ui,...,Um, as a 
polynomial function q, in some pre-processed parameter dependent values 
ui{ui, ...,Um), ...jOJi^ui, ...,Um) such that, for the evaluation of p in the input 
y, no branchings are needed anymore. 

In the language of constraint databases, this pre-processing of 
) into the functions a;i(iti, ^^(wi, m^), can 
be modeled as a query with input schema S = (G) and output schema 
S = {Qi, ...,Q£), with ar{Qj) = m {j = 1, ...,£). We remark that, as soon as 
ui, ...,Um and i are fixed, also the database schema S becomes fixed. 

In the geometric elimination literature, algorithms are (implicitly or ex- 
plicitly) required to be robust (see [7j for the definition of the notion of 
robustness). In the formalism of constraint databases, robustness can be 
specified by requiring that the query, that transforms g{ui, Um, xi, . . . , Xn) 
into ui{ui, ...,Um), ■■■,0Je{ui, ...,Um), is a geometric query. This means 
that for any two parameter instances {ai,...,am,) and {a[,...,a'^) for 
which g{ai, ...,am,Xi, . . . ,Xn) and g{a[, a'^, xi, . . . , Xn) determine the 
same polynomials in xi,...,Xn, also a'i(Q;i, a^) = '^i(o'i; <^m); •••^ 
coi^ai, ...,am) = uje{a[, a'^) holds. Therefore we may conclude that the 
query that transforms the input database instance (g) into the output data- 
base instance [p) is also a geometric query. On the other hand, we observe 
that u{ui, Urn), '^^('Ui, ...Um) may be interpreted as parameters of a fixed 
division-free circuit with input y, which is represented by the polynomial q 
and which is independent of the input g. This circuit evaluates the polyno- 
mial p{Ui, ...,Um,y)- 

The task of transforming the elimination problem (*a), given by the poly- 
nomial g{ui, ...,Um,Xi, ...,Xn), iuto a collection of intermediate polynomi- 
als uJi{ui, ...,Um), ...,Ui{ui, ...,Um) that satisfy a number of restrictions, can 
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therefore be specified using the constraint database formahsm, as we have 
just illustrated. 

We conclude this section by mentioning that we may ask how long the 
vector {ui{ui, ...,Um), ■.■,u}e{ui, ...,Um)) has to be in order to be able to ex- 
press p{ui, ...,Um,y) by a fixed division-free arithmetic circuit with input y, 
which is independent of the polynomial g. As the main result of this paper, 
it turns out that the notion of geometric query implies, that the number i, 
necessarily has to be exponentially big in the number of quantified variables, 
i.e., i > 2". 

4 Representation of databases and algorith- 
mic model for query evaluation 

4.1 Representation of databases 

The function symbols appearing in a database schema will be interpreted by 
polynomial functions that are well-defined in any commutative ring. Typi- 
cally, the data structure to implement these functions are division-free arith- 
metic circuits that allow the evaluation of the function in an arbitrary com- 
mutative ring. 

The relation symbols appearing in a database schema will be interpreted 
by relations that may be evaluated in an arbitrary commutative ring. The 
relations are implemented by arithmetic boolean circuits. 

Section [2] 

These representations may be considered as implementations of a certain 
datatype which allows in case of function symbols the evaluation of the func- 
tion in an arbitrary commutative ring and in case of relation symbols the 
evaluation of membership 

This paper uses notions from algebraic geometry and commutative alge- 
bra. These notions can be found in standard textbooks, such as [22] and we 
refer to Appendix A for an overview of them. 
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4.2 Algorithmic models and complexity measures. 

4.2.1 Data structures: essentially division- free arithmetic circuits 

The algorithmic problems considered in this paper will depend on continu- 
ous parameters and therefore the corresponding input data structures have 
to contain entries for these parameters. We call them problem or input pa- 
rameters. Once such a parametric problem is given, the specialization of 
the parameters representing input objects are called (admissible) problem or 
input instances and these may in principle be algebraically dependent. An al- 
gorithm solving the given problem operates on the corresponding input data 
structure and produces for each admissible input instance an output instance 
which belongs to a previously chosen output data structure. 

The procedures or algorithms considered in this paper operate with es- 
sentially division-free arithmetic circuits as basic data structures for the rep- 
resentation of inputs and outputs. An essentially division-free arithmetic 
circuit (or straight-line program) is an algorithmic device that can be rep- 
resented by a labeled directed acyclic graph {dag) as follows: the circuit 
depends on certain input nodes, labeled by indctcrminatcs over the field 
of the rationals Q. These indeterminates are thought to be subdivided in 
two disjoints sets, representing the parameters and the variables of the given 
circuit. For the sake of definiteness, let C/i, . . . , C/^ be the parameters and 
Yi, . . . ,Yf the variables of the circuit. Let K := Q{Ui, . . . , Ur) be the field 
of quotients of polynomials in Q[f/i, . . . , U,]. Wc call K the parameter fi,eld 
of the circuit. The circuit nodes of indegree zero which are not inputs are 
labeled by elements of Q, which are called the scalar s of the circuit (here 
"indegree" means the number of incoming edges of the corresponding node) . 
Internal nodes are labeled by arithmetic operations (addition, subtraction, 
multiplication and division). We require that the internal nodes of the circuit 
represent polynomials in the variables Yi, . . . , 1^. We call these polynomials 
the intermediate results of the given circuit. The coefficients of these polyno- 
mials belong to the parameter field K. In order to achieve this requirement, 
we allow in an essentially division-free circuit only divisions which involve 
elements of K. Thus essentially division-free circuits do not contain divi- 
sions involving intermediate results which depend on the variables Yi, . . . , 1^. 
A circuit which contains only divisions by non-zero elements of Q is called 
totally division- free. 

Finally we suppose that the given circuit contains one or more nodes 
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which are labeled as output nodes. The results of these nodes are called 
outputs of the circuit. Output nodes may occur labeled additionally by sign 
marks of the form "= 0" or "7^ 0" or may remain unlabeled. Thus the 
given circuit represents by means of the output nodes which are labeled by 
sign marks a system of parametric polynomial equations and inequations. 
This system determines in its turn for each admissible parameter instance a 
locally closed set (i.e., an embedded affine variety) with respect to the Zariski 
topology of the afiine space A* of variable instances. The output nodes of the 
given circuit which remain unlabeled by sign marks represent a parametric 
polynomial application (in fact a morphism of algebraic varieties) which maps 
for each admissible parameter instance the corresponding locally closed set 
into a suitable affine space. We shall interpret the system of polynomial 
equations and inequations represented by the circuit as a parametric family of 
systems in the variables of the circuit. The corresponding varieties constitute 
a parametric family of varieties. The same point of view is applied to the 
morphism determined by the unlabeled output nodes of the circuit. We shall 
consider this morphism as a parametric family of morphisms. 

4.2.2 Complexity models 

To a given essentially division-free arithmetic circuit we may associate dif- 
ferent complexity measures and models. In this paper we shall be mainly 
concerned with sequential computing time^ measured by the size of the cir- 
cuit. Occasionally we will also refer to parallel time., measured by the depth 
of the circuit. In our main complexity model is the total one, where we take 
into account all arithmetic operations (additions, subtractions, multiplica- 
tions and possibly occurring divisions) at unit costs. For purely technical 
reasons we shall also consider two non-scalar complexity models, one over 
the ground field Q and the other one over the parameter field K. In the non- 
scalar complexity model over K we count only the essential multiplications 
(i.e. multiplications between intermediate results which actually involve vari- 
ables and not only parameters). This means that linear operations (i.e. 
additions and multiplications by arbitrary elements of K) are cost free. Sim- 
ilarly, Q-linear operations are not counted in the non-scalar model over Q. 

Let ^1, . . . , 6*^ be the elements of the parameter field K computed by the 
given circuit. Since this circuit is essentially division-free we conclude that 
its outputs belong to Q[^i, . . . , ^m][^i) ■ ■ ■ i Let L be the non-scalar size 
(over K) of the given circuit and suppose that the circuit contains q output 
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nodes. Then the circuit may be rearranged (without affecting its non-scalar 
complexity nor its outputs) in such a way that the condition 

m = L'^ + {2t-l)L + q{L + t + l) (1) 

is satisfied. In the sequel we shall always assume that 

we have already performed this rearrangement. Let 
Y := 9 := and let /i,...,/, G Q[e][Y] be the 

outputs of the given circuit. Let Zi, . . . , Zm be new indeterminates and write 
Z := (Zi, . . . , Zra)- Then there exist polynomials Fi, . . . , Fg G Q[Z, Y] such 
that /i = y), ...,/, = Fgie, Y) holds. Let us write / := (A, . . . , /,) 
and F := (Fi, . . . , Fg). Consider the object class 

O := {F{C,Y) : C e A™} 

which we think represented by the data structure V := A™ by means of the 
obvious encoding which maps each code ( GVto the object F{(, Y) G C[y]^. 

For the moment, let us consider as input data structure the Zariski open 
subset U C A*" where the rational map 9 = {9i, . . . ,9^) is defined. Then 
the given essentially division-free arithmetic circuit represents an algorithm 
which computes for each input code w G W an output code 9{u) representing 
the output object f{u,Y) = F(^6(u),Y). This algorithm is in the above 
sense essentially division-free. From identity ([1]) we deduce that the size m 
of the data structure V is closely related to the non-scalar size L of the given 
circuit. In particular we have the estimate 

Vm-{t + q)<L. (2) 

Later we shall meet specific situations where we are able to deduce 
from a previous (mathematical) knowledge of the mathematical object 
/ = (/i, . . . ,/q) a lower bound for the size of the output data structure 
of any essentially division-free algorithm which computes for an arbitrary 
input code u eU the object f{u,Y). Of course, in such situations we ob- 
tain by means of ([2]) a lower bound for the non-scalar size (over K) of any 
essentially division-free arithmetic circuit which solves the same task. In 
particular we obtain lower bounds for the total size and for the non-scalar 
size over Q of all such arithmetic circuits. 

4.2.3 Elimination problems 

Given an essentially division- free arithmetic circuit as input, an elimination 
problem consists in the task of finding an essentially division-free output 
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circuit which describes the Zariski closure of the image of the morphism 
determined by the input circuit. The output circuit and the corresponding 
algebraic variety are also called a solution of the given elimination problem. 
We say that a given parameter point fixes an instance of the elimination 
problem under consideration. In this sense a problem instance is described 
by an input and an output (or solution) instance. 

In this paper we restrict our attention to input circuits which are totally 
division-free and contain only output nodes labelled by "=0" and unlabelled 
output nodes. Mostly our output circuits will also be totally division- 
free and will contain only one output node, labelled by the mark "=0". 
This output node will always represent a canonical elimination polynomial 
associated to the ehmination problem under consideration (see Section 14.31 
for more details). 

In case that our output circuit contains divisions (depending only on 
parameters but not on variables), we require to be able to perform these 
divisions for any problem instance. In order to make this requirement sound, 
we admit in our algorithmic model certain limit processes in the spirit of de 
I'Hopital's rule (below we shall modelise these limit processes algebraically, 
in terms of places and valuations). The restriction we impose on the pos- 
sible divisions in an output circuit represents a first fundamental geometric 
uniformity requirement for our algorithmic model. 

4.2.4 Parametric elimination procedures 

An algorithm which solves a given elimination problem may be considered 
as a (geometric) elimination procedure. However this simple minded notion 
is too restrictive for our purpose of showing lower complexity bounds for 
elimination problems. It is thinkable that there exists for every individual 
elimination problem an efficient ad hoc algorithm, but that there is no uni- 
versal way to find and to represent all these ad hoc procedures. Therefore, a 
geometric elimination procedure in the sense of this paper will satisfy certain 
uniformity and universality requirements which we are going to explain now. 

We modelise our elimination procedures by families of arithmetic 
networks (also called arithmetic-boolean circuits) which solve entire 
classes of elimination problems of arbitrary input size. In this sense we shall 
require the universality of our geometric elimination procedures. Moreover, 
we require that our elimination procedures should be essentially division-free. 

In a universal geometric elimination procedure, branchings and divisions 
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by intermediate results (that involve only parameters, but not variables) 
cannot be avoided. From our elimination procedures we shall require to 
be parsimonious with respect to branchings (and divisions). In particular 
we shall require that our elimination procedures do not introduce branchings 
and divisions for the solution of a given elimination problem when traditional 
algorithms do not demand this (an example of such a situation is given by 
the flat families of elimination problems we are going to consider below in 
Section 14. 3p . This restriction represents a second fundamental uniformity 
requirement for our algorithmic model. 

We call a universal elimination procedure parametric if it satisfies our first 
and second uniformity requirement, i.e., if the procedure does not contain 
branchings which otherwise could be avoided and if all possibly occurring 
divisions can be performed on all problem instances, in the way we have 
explained before. In this paper we shall only consider parametric elimination 
procedures. 

We call a parametric elimination procedure geometrically robust if it pro- 
duces for any input instance an output circuit which depends only on the 
mathematical objects "input equation system" and "input morphism" but 
not on their circuit representation. We shall apply this notion only to elimi- 
nation problems given by (geometrically or scheme-theoretically) fiat families 
of algebraic varieties. This means informally that a parametric elimination 
procedure is geometrically robust if it produces for fiat families of problem 
instances "continuous" solutions. 

Of course, our notion of geometric robustness depends on the (geometric 
or scheme-theoretical) context, i.e. it is not the same for schemes or varieties. 
In Section H73] we shall explain our idea of geometric robustness in the typical 
situation of fiat families of algebraic varieties given by reduced complete 
intersections. 

Traditionally, the size of a system of polynomial equations (and inequa- 
tions) is measured in purely extrinsic, syntactic terms (e.g. number of pa- 
rameters and variables, degree of the input polynomials, size and depth of 
the input circuit etc). However, there exists a new generation of symbolic 
and numeric algorithms which take also into account intrinsic, semantic (e.g. 
geometric or arithmetic) invariants of the input equation system in order to 
measure the complexity of elimination procedure under consideration more 
accurately. 

In this paper we shall turn back to the traditional point of view. In 
[7] it was shown that, under certain universality and uniformity restrictions. 
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no parametric elimination procedure which includes efficient computation of 
Zariski closures and of generically squarefree parametric greatest common 
divisors for circuit represented algebraic families of polynomials, is able to 
solve an arbitrary elimination problem in polynomial (sequential) time, if 
time is measured in terms of circuit size and input length is measured in 
syntactical terms only. 

4.3 Flat families of elimination problems 
4.3.1 Definitions and preliminaries 

Let, as before, let Ui, . . . ,Ur, Xi, . . . , X„, Y be indeterminates over Q. In the 
sequel we shall consider Xi, . . . , Xn and Y as variables and Ui, . . . ,Ur as pa- 
rameters. Let Gi, . . . ,Gn and F be polynomials belonging to the Q-algebra 
Q[[/i, . . . , Ur, Xi, . . . , Xn]- Suppose that the polynomials Gi, . . . , Gn form a 
regular sequence in Q[Ui, . . . ,Ur, Xi, . . . , Xn] defining thus an equidimen- 
sional subvariety V := {Gi = 0, . . . ,Gn = 0} of the (r + n)-dimensional 
affine space A'' x A" over the field C. The algebraic variety V has di- 
mension r. Let S be the (geometric) degree of V (this degree does not 
take into account multiplicities or components at infinity). Suppose fur- 
thermore that the morphism of affine varieties tt : V — > A*", induced by the 
canonical projection of A*" x A" onto A*", is finite and generically unramified 
(this implies that tt is flat and that the ideal generated by Gi, . . . ,Gn in 
Q[Ui, . . . ,Ur.,Xi, . . . , Xn] is radical). Let tt : V — > A''+^ be the morphism 
defined by tt{z) := (tt^z), F{z)) for any point z of the variety V. The im- 
age of TT is a hypersurface of A''+^ whose minimal equation is a polynomial 
of k[Ui, . . . , Ur,Y] which we denote by P. Let us write degP for the total 
degree of the polynomial P and degy P for its partial degree in the variable 
Y. Observe that P is monic in Y and that deg P < 5 deg F holds. Further- 
more, for a Zariski dense set of points it of A*", we have that degy P is the 
cardinality of the image of the restriction of F to the finite set 'k~^{u). The 
polynomial P{Ui, . . . , Ur, F) vanishes on the variety V. 

Let us consider an arbitrary point u — {ui, . . . ,Ur) of A''. For ar- 
bitrary polynomials A e Q[Ui, . . .,Ur, Xi, . . ., Xn] and B e Q[Ui, . . . ,Ur,Y] 
we denote by A^'^^ and B^'^^ the polynomials A{ui, . . . ,Ur, Xi, . . . , Xn) 
and B{ui, . . . ,Ur,Y) which belong to Q('Ui, • • • , Ur)[Xi, .... X„] and 
Cl{ui, . . . ,Ur)[Y] respectively. Similarly we denote for an arbitrary poly- 
nomial C e Q[Ui, . . . ,Ur] by C^"^ the value C{ui, . . . ,Ur) which belongs to 
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the field Q(ui, . . . , Ur). The polynomials \ . . . , Gn define a zero dimen- 
sional subvariety y(^) := {gS"^ = 0, . . . , Ci"^ = 0} = n-\u) of the affine 
space A". The degree (cardinality) of V^"^^ is bounded by 6. Denote by 
TT : V^(") — y the morphisms induced by the polynomial F^^^ on the va- 
riety Observe that the polynomial P^^^ vanishes on the (finite) image of 
the morphism tt^"^. Observe also that the polynomial P*^"-* is not necessarily 
the minimal equation of the image of tt^"^. 

We call the equation system Gi = 0, . . . , = and the polynomial F a 
flat family of elimination problems depending on the parameters Ui, . . . ,Ur 
and we call P the associated elimination polynomial. An element u G A'' 
is considered as a parameter point which determines a particular problem 
instance (see Section 14^2]) . The equation system Gi = 0, . . . , G„ = together 
with the polynomial F is called the general instance of the given fiat family of 
elimination problems and the elimination polynomial P is called the general 
solution of this fiat family. 

The problem instance determined by the parameter point u G A'' is given 
by the equations G["^ = 0,...,G-n^ = and the polynomial The 
polynomial P^"-* is called a solution of this particular problem instance. We 
call two parameter points u, u' G A'' equivalent (in symbols: u ~ u') if 
^{u) ^ ^{n')^ . . . ^ ci"^ = G^n^ and = F^"') holds. Observe that u u' 
implies F^") = We call polynomials A G Q[[/i, . . . , [/^, Xi, . . . 

-B G Q[[/i, . . . , f/r, Y] and G G Q[t/i, . . . , f/r] invariant (with respect to ~) if 
for any two parameter points u,u' of with u ^ u' the respective identities 
_ ^K)^ = ^K) and G(") = G^''') hold. 

4.3.2 Arithmetic circuits 

An arithmetic circuit in Q[[/i, . . . , [4, F] with scalars in Q[[/i, . . . , [4] is a 
totally division-free arithmetic circuit in Q[t/i, . . . , Ur, Y], say /3, modehsed in 
the following way: (5 is given by a directed acyclic graph whose internal nodes 
are labelled as before by arithmetic operations. There is only one input node 
of /3, labelled by the variable Y . The other nodes of indegree zero the circuit 
(5 may contain, are labelled by arbitrary elements of Q[t/i, . . . , f/r]- These 
elements are considered as the scalars of (3. We call such an arithmetic circuit 
13 invariant (with respect to the equivalence relation ~) if all its scalars 
are invariant polynomials of Q[Ui, . . . ,Ur]. Considering instead of Y the 
variables Xi, . . . , X„ as inputs, one may analogously define the notion of an 
arithmetic circuit in Q[[/i, . . . , [4, Xi, . . . , X„] with scalars in Cl[Ui, . . . ,Ur] 
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and the meaning of its invariance. However, typically we shall limit ourselves 
to circuits in Q[Ui, . . . , Ur, Y] with scalars in Q[Ui, . . . ,Ur]- 

4.3.3 Geometrically robust parametric elimination problems 

We are now ready to characterise in the given situation what we mean by 
a geometrically robust parametric elimination procedure. Suppose that the 
polynomials Gi, . . . , Gn and F are given by a totally division-free arithmetic 
circuit /3 in Q[C/i, . . . , C/^., Xi, . . . , X^]. A geometrically robust parametric 
elimination procedure accepts the circuit /? as input and produces as out- 
put an invariant circuit F in Q[L^i, . . . , Ur, Y] with scalars in Q[Ui, . . . , Ur], 
such that r represents the polynomial P. Observe that in our definition of 
geometric robustness we did not require that /3 is an invariant circuit be- 
cause this would be too restrictive for the modelling of concrete situations in 
computational elimination theory. 

The invariance property required for the output circuit F means the fol- 
lowing: let u = {ui, . . . ,Ur) be a parameter point of and let F*^"^ be the 
arithmetic circuit in Q(«i, . . . , obtained from the circuit F evaluating 

in the point u the elements of Q\Ui, . . . ,Ur] which occur as scalars of F. Then 
the invariance of F means that the circuit F^") depends only on the particular 
problem instance determined by the parameter point u but not on u itself. 
Said otherwise, a geometrically robust elimination procedure produces the 
solution of a particular problem instance in a way which is independent of 
the possibly different representations of the given problem instance. 

By definition, a geometrically robust parametric elimination procedure 
produces always the general solution of the flat family of elimination prob- 
lems under consideration. This means that for flat families, geometrically 
robust parametric elimination procedures do not introduce branchings in the 
output circuits. It turns out that the following meta-statement becomes 
true: within the standard philosophy of commutative algebra, none of the 
known (exponential time) parametric elimination procedures can be improved 
to a polynomial time algorithm. For this purpose it is important to remark 
that the known parametric elimination procedures (which are without ex- 
ception based on linear algebra as well as on comprehensive Grobner basis 
techniques) are all geometrically robust. 

The invariance property of these procedures is easily verified in the situa- 
tion of a fiat family of elimination problems. One has only to observe that all 
known elimination procedures accept the input polynomials Gi, . . . , G„ and 
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F in their dense or sparse coefficient representation or as evaluation black 
box with respect to the variables Xi, . . . , Xn- 

Finally let us observe that robust elimination procedures can be specified 
as geometric queries. 

5 A lower complexity bound for evaluating 
geometric queries 

As main result of this paper, we obtain that for geometric FO - queries the 
size of the output data schema may become necessarily exponentially big in 
the number of quantified variables occurring in the query. 

The proof of this fact goes along the lines of Theorem below. 

6 Sample point queries 

6.1 Sample point queries and generalized sample point 
queries 

6.1.1 The example of rationally pcirameterized families of polyno- 
mial functions. 

A particular instance of interest is the case that the semi- algebraic set A is 
contained in and represents a rational family of polynomial functions 

from R" to R. To be more precise, let tt : R"^+"+i R™ be the canonical 
projection of any point of R"^+"+i on its first m coordinates. Suppose that A 
is non-empty and that for any u — . . . , Um) £ the semi- algebraic set 
({■u} X R"+^) n A is the graph of an n-variate polynomial /„ G R[a;i, . . . , a;„]. 
It is a natural extension of our previously introduced sample-point query to 
ask for a procedure which enables us for each u G '^{A) and each x G R" 
to compute the value of fu{x). The output of such a procedure may be a 
purely existential prenex first-order formula in the free variables Ui, . . . ,u.m 
and which represents for each u G 7t{A) a division- free arithmetic 

circuit which evaluates the polynomial /„ (observe that there exists a uni- 
form degree bound for all these polynomials). One easily verifies that all 
our requirements on the semi-algebraic set A, except that of the polynomial 
character of the function represented by the graph {{u} x R"+^) flA, are first- 
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order definable over the reals. Nevertheless, over the complex numbers, when 
A is a constructible (i.e., a first-order definable subset of C'""'""'"'"^), all these 
requirements are first-order expressible. This leads us to a new type of com- 
putable queries which return on input a semi-algebraic or constructible set 
A as above and an element u G it {A), a first-order formula which represents 
a division-free arithmetic circuit evaluating the polynomial /„. Uniformity 
of query evaluation with respect to u is expressed by the requirement that 
the terms contained in this formula have to depend rationally on u. 

In the following, we shall refer to this type of queries as extended sample 
point queries. We shall refer to ui, ... ,Um 8.s the parameters and to xi, . . . ,Xn 
as the variables of the query. 

Suppose now that A is a constructible subset of C"^+"-+^ with irreducible 
Zariski-closure. Let V be the Zariski-closure of tt{A). Then V is an irre- 
ducible affine subvariety of C™. We denote the function field of V by K. It 
is not difficult to see that for generically chosen parameter points u G 7^{A), 
the extended sample point query associated to A can be realized by a greatest 
common divisor computation in the polynomial ring K[xi, . . . , Xn]. 

6.1.2 Variables versus parameters. 

The previous example is motivated by spatial data that come from physical 
observation and are only known with uncertainty. Another motivation comes 
from parametric optimization. Optimization problems described in can also 
be studied in parametric form, i.e., in the case where the linear inequalities 
and the target function contain coefficients that depend on a time parameter 
[5j and [H] . In this case, an optimum is not an arbitrary set of sample points 
but an analytic (or at least continuous) function which depends on a time 
parameter. 

We are now going to explain why we distinguished between the parameters 
Ui, . . . ,Um and the variables Xi, . . . , 2;„ in our discussion of rationally param- 
eterized families of polynomial functions. In the example above, let $(wi, 
. . . ,Um', xi, . . . ,Xm,y) be a quantifier-free formula which defines the semi- 
algebraic or constructible set A. Let us suppose that $ contains a subformula 
\1/(mi, . . . , Urn) which expresses an internal algebraic dependency between the 
parameters Ui, . . . , Um- With respect to the variables Xi, . . . , x„ there is no 
such subformula contained in $. For the sake of simplicity, we shall suppose 
that \1/ defines the set 'k{A) and that there exists a formula VL{ui, . . . , Um', a^i, 
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. . . , Xn, y) such that $(iti, . . . , Um] Xi, . . . , Xjn-i v) Can be written as 

...,Um) a ^{Ui, ...,Um;Xi,.. .,Xm,y)- 

Below we shall meet natural examples of parameterized algebraic families 
of polynomial functions where the formula ^' becomes of uncontrolled size 
and is of few interest, whereas the formula Q becomes the relevant part of the 
output information of a suitable elimination algorithm. This situation occurs 
for instance when the points u of R"* satisfying the formula ^ are given in 
parametric form (i.e., when they are image points of some polynomial or 
rational map coming from some affine source space). In this case, we are 
only interested in the subformula Q, since points u e R"* satisfying ^ can 
easily be produced in sufficient quantity. In subsequent queries, Xi, . . . ,x„ 
may appear as bounded variables, whereas the parameters Ui,...,Um are 
not supposed to be subject to quantification. The example of the Ui, . . . , Um 
expressing uncertainty in physical spatial data illustrates this. These different 
roles motivate us to distinguish between ui,.. . ,Um and and to call 

them parameters and variables, respectively. 

6.1.3 The branching-parsimonious algorithmic model. 

In the model, that we are going to use in the sequel, parameters and variables 
receive a different treatment. (Free) variables may be specialized into arbi- 
trary real (or complex) values, whereas the specialization of parameters may 
be subject to certain restrictions. In the above example of a rationally param- 
eterized family of polynomial functions, the n-tuple of variables (xi, . . . , x„) 
may be speciahzed into any point of the affine space R" (or C"), whereas the 
m-tuple of parameters (iti, . . . , Um) may only be specialized into points sat- 
isfying \E'('Ui, . . . ,«„), i.e., into points belonging to 7r(A). Once the ri-tuple 
of variables (xi, . . . , Xn) is specialized into a point of the corresponding affine 
space, this point cannot be modified anymore. However, we allow infinitesi- 
mal modifications of a given specialization of the m-tuple of parameters (iti, 
. . . , Ujn) within the domain of definition determined by the formula ^ . In 
the branching-parsimonious model, we require that an arithmetic boolean 
circuit which represents the semi-algebraic or constructible set A does not 
contain divisions which involve the variables xi, . . . , a;„. Similarly, for a given 
point u e 77(^4), we require that the arithmetic circuit representing the poly- 
nomial fu is division-free. However, divisions by algebraic expressions in 
the parameters iti,...,^^ are sometimes unavoidable (e.g., in the case of 
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parametric greatest common divisor computations; see [Hin]). Therefore, we 
allow certain limited divisions by algebraic expressions which depend only 
on the parameters ui, . . . , Um- More precisely, we allow that the arithmetic 
boolean circuits representing the set A or the output of the corresponding 
extended sample point query computes certain, but not arbitrary, rational 
functions in the parameters Ui, . . . , Um, called scalars of the circuit. However, 
we do not allow the division of a positive- degree polynomial in the variables 
Xi, . . . , by a non-constant scalar. In the above sense, we require for our 
branching-parsimonious algorithmic model that arithmetic boolean circuits 
are essentially division-free with respect to variables (see [H [9] for a precise 
definition). 

6.1.4 Branching- free output representations of extended sample 
point queries. 

Since we allow certain infinitesimal modifications of the parameters lii, . . . , 
Um within their domain of definition, we sometimes may replace divisions 
(and corresponding branchings) by limit processes in the spirit of L'Hopital's 
rule. It is possible to mimic algebraically this kind of limit process by places 
(see |19j for the notion of place and [3 [H] for motivations of this idea). 

Branchings corresponding to divisions can trivially be avoided by restrict- 
ing input data. Therefore a meaningful notion of branching-parsimonious (or 
branching-free) algorithm requires the consideration of Zariski-closures of in- 
put data sets. This may partially explain the rather technical assumptions 
and tools in the following ad hoc definition of a branching-free representation 
of the output of an extended sample point query. 

Suppose now that in the example above A is a constructible subset of 
Qm+n+i ^j^]-^ irreducible Zariski- closure B. Let V be the Zariski-closure of 
7r(yl) in C™. Then V is an irreducible afiine variety whose function field 
we denote by K. Moreover, the irreducible affine variety B is birationally 
equivalent to y x C". Suppose furthermore that 7r(5) = V holds and that B 
represents a rationally parameterized family of polynomial functions which 
extends the family represented by A. Then we say that the extended sample 
point query associated with A admits a branching-free output representation 
if there exists an essentially division-free, single-output arithmetic circuit /3 
with inputs and scalars 9i, ... ,9s G K satisfying the following 

conditions: 

(i) for any point u G tt{A) where the rational functions 9i,...,6s are 
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defined, the division-free aritlimetic circuit, obtained from /3 by spe- 
cializing tlie scalars 6i, ... ,6s into the complex values 9i{u), . . . , 9s{u), 
evaluates the polynomial /„; 

(ii) for any point u & V and any place ip : K — )■ CU{oo} whose valu- 
ation ring extends the local ring of the affine variety V at the point 
u, the values v^(^^i), . . . , '^{ds) are finite and uniquely determined by u 
(therefore we shall write 6i{u) := ip{9i), . . . ,9s{u) := tpifis))- 

Let an arithmetic circuit /3 be given as above. Then we call (3 a branching- 
free representation of the output of the extended sample point query associ- 
ated to A. 

Observe that the output of the circuit (3 represents a polynomial belong- 
ing to K[xi, . . . whose coefficients satisfy condition (ii). Moreover, the 
arithmetic circuit /3 constitutes a division-free representation of the extended 
sample point query associated to the Zariski-closure B of A. Finally, let us 
remark that for any w e y, x G C" and y G C, the point {u,x,y) belongs 
to B if and only if the circuit ^u, obtained from (3 by replacing the scalars 
9i, ... ,9s by the complex numbers 9i{u), . . . , 9s{u), computes on input x the 
output y. 

We require that a branching-parsimonious query evaluation algorithm 
produces a branching-free output representation of the given extended sam- 
ple point query if the query admits such a representation. 

Let us also observe that extended sample point queries appear in a natural 
way if we apply the constraint database concept to data processing in the 
context of approximation theory and functional analysis. 

6.2 A lower complexity bound for evaluating sample 
point queries 

In this section, we restrict our attention to constraint databases defined in the 
language F0(-|-, x, 0, 1, =) over the complex numbers. We shall consider two 
ternary relational predicates, namely S{vi,V2, w) and P{vi, V2,w). Our query 
language will therefore be F0(-|-, x, 0, 1, =, 5', P). Let L, n be given natural 
numbers and let r :— {L + n + 1)^. For any polynomial / e C[xi, . . . ,Xn], 
we denote by L{f) the minimal non-scalar size of all division-free arithmetic 
circuits with inputs xi, . . . ,Xn and scalars from C which evaluate the poly- 



35 



nomial /. Let 

WL,n--={f eC[x,,...,Xr.]\L{f)<L}. 

One sees easily that all polynomials contained in WL,n have degree at most 
2^ and that WL,n forms a Q-definable object class which has a Q-definable 
holomorphic encoding by the continuous data structure C. Observe that 
Zariski-closure Wi^n of WL,n is a Q-definable, absolutely irreducible algebraic 
variety consisting of the polynomials of C[xi, . . . , Xn] which have approximate 
complexity at most L. Moreover, the affine variety Wi^n forms a cone in its 
ambient space (i.e., for any A G C we have XWL,n ^ WL,n)- For details on 
complexity and data structure models we refer to [H [7] . 

Let zi, . . . , Zr and y be new variables. Choose now a directed acyclic graph 
Vi^n representing a generic, division-free arithmetic circuit with input nodes 
Xi, . . . ,Xn, output node y and scalar nodes zi, . . . , Zr such that any polyno- 
mial of W^^n may be evaluated by the division-free arithmetic circuit obtained 
from by a suitable specialization of the parameters zi, . . . ,Zr into com- 
plex values. Without loss of generality, we may assume that the number of 
internal nodes of Vi^n is of order 0{{L + n)^). Translating the structure 
of the directed acyclic graph Vi^n into first-order logic one infers easily a 
formula \I/l^„(S', P, zi, . . . , Zr,Xi, . . . , Xn,y) in the free variables zi, . . . , Zr,Xi, 
. . . , Xn, y of the query language F0(+, x , 0, 1, =, S, P) such that L,n{S, P) 
satisfies the following conditions. 

(i) \&L,„(S', P) is prenex, purely existential and of length 0((L + n)^); 

(ii) interpreting the predicates S and P in ^ L,n as the graphs of the addi- 
tion and the multiplication of complex numbers, and specializing the 
variables zi, . . . ,Zr into the complex numbers Ci, . . . , (^r, the formula 
^L^„(S', P, Ci, • • • , Cr> 3^1, ■ ■ ■ , describes the graph of the polyno- 
mial of C[xi, . . . , Xn] computed by the arithmetic circuit, obtained from 
T^L,n by specializing the scalars zi, . . . ,Zr into Ci, • • • , Cr- 

Let m := 4(L -|- n)^ + 2. From [7], Corollary 3 (see also [9i Lemma 4]) we 
deduce that there exists an identification sequence 71, ... ,7^ G Q" for the 
object class VFl,„. Let Aj, „(S', P) be a closed F0(+, x , 0, 1, =, 5", P)-formula 
saying that S and P are the graphs of two binary operations which map 

into C, that 71, . . . ,7„ is an identification sequence for the object class 
of applications from C" to C, defined by the F0(+, x , 0, 1, =, S", F)-formula 
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^L,niS, P) and that this object class is not empty. Without loss of generality, 
we may assume that A^, „(S', P) has length 0{{L + n)^) and is prenex with 
a fixed number of quantifier alternations (which is independent of L and n). 

We consider now the F0(+, x, 0, 1, =, S*, P)-formulas ^L,n{S, P,ui, . . . , 
Ur,xi, . . . , Xn, y) defined by 

{3Zi) ■ ■ ■ {3Zr){^L,niS, P, Zi, . . . , Zr, Xi, . . . , Xn, y)A 

/\ ^L,n('S', P,Zi,..., Zr, 7fe, Mjfc)) A Al,„(S', P) 
l<k<m 

and flL,n{S, P,Ui, . . . , u^) defined by 

{3zi) ■ ■ ■ {3Zr){ /\ *L,„(S,P,Zi,...,Z^,7fe,life)) AAl,„(5,P). 
l<fe<m 

Without loss of generality, we may assume that ^L^n{S, P) and Q,L^n{S, P) 
are prenex formulas of length 0{{L + n)^) having a fixed number of quanti- 
fier alternations and containing the free variables Mi, . . . , Um', Xi, . . . , Xn, y and 
Ml, . . . ,Um, respectively. Let tt : C™'"^""'"^ — ?■ C™" be the canonical projection 
which maps each point of C'""'"""'"^ on its first m coordinates and let D be a 
constraint database over the schema {S, P) over the complex numbers. Sup- 
pose that D satisfies the formula Ai^„(5', P). With respect to the database 
D, the formula ^L,n{S,P, y) defines a non-empty con- 

structible subset Ai^niD) of and the formula QL,n{S, P,ui, . . . , Ur) 

defines the set 7r{AL^n{D))- Moreover, for any u G TT{AL,n{D)), the for- 
mula ^L^n{S, P,u,Xi, . . . ,Xn,y) describes the graph of a n-variate polyno- 
mial map. Therefore, it makes sense to consider, for any natural numbers n 
and L, the generalized sample point query associated to the formula ^L^n{S, 
P,ui, . . . ,Un,Xi, . . . ,Xn,y)- Supposc uow that there is given a branching- 
parsimonious procedure V which evaluates this family of extended sample 
point queries. We are now going to analyze the complexity behaviour of 
V for this query on the particular input database D, where S and P are 
interpreted as the graphs of the sum and the product of complex numbers. 

We are now able to state and to prove the main complexity result of this 
paper. 

Theorem 1 Let notations and assumptions be as before. Then the 
branching-parsimonious procedure V requires sequential time 2^^'^^ in order 
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to evaluate on input the database D the extended sample point query associ- 
ated to the size 0{n?) first-order formula $„^„(S', P). In particular, extended 
sample point queries associated to first-order formulas with a fixed number of 
quantifier alternations cannot be evaluated by branching-parsimonious proce- 
dures in polynomial time. 

Proof. The arguments we are now going to use follow the general lines of 
the proofs of [9, Theorem 5] and [3 Theorem 4]. 

For the moment let us fix the integer parameters L and n. 

Observe that the closed formula Al „(S', F) is valid on the database D. 
Therefore the constructible set A^^n '■= ^L,n(-D) is nonempty. 

Let Bi n and Vi^n be the Zariski-closures of A^^n and Tr^A^ n) in C"^+"-+^ 
and C™, respectively. 

Let XL,n ■■= WL,n X C" ^ 0™+"+^ and ^L,n ■■= WL,n C"^ be the 
morphisms of Q-definable affine varieties defined for / G Wi^n and x G C" by 

^L,n{f^^) ■= ifill),---Jhm),xJ{x)) and /iL,n(/) := (/(7l),---,/(7m))- 

From the syntactic form of $l^„(6', P) and Q^^^iS, P) one infers immediately 
that 

XL,n{WL,n X C") = A^,, 

and 

fJ^L,n(WL,n) = Tr{AL,n) 

holds. 

Therefore B^ n and V^^n are Q-definable absolutely irreducible affine va- 
rieties and we may consider „ and as dominant morphism mapping 
W^L,n X C" into B^ n and W^l,„ into VL,n- Observe that WL,n and V^^n form 
closed cones in their respective ambient spaces. Since 71, . . . , 7^ were chosen 
as an identification sequence for the object class W^^n, we may conclude that 
fJ'L,n '■ W L,n — ^ Vl^u is an injective dominant morphism of closed affine cones, 
which is homogeneous of degree one. 

Therefore iiL,n is a finite, bijective and birational morphism of affine va- 
rieties (see, e.g., [221 L5.3 Theorem 8 and proof of Theorem 7], [9l Lemma 4] 
or lU Lemma 5]). 

Let TT : C™'"'""^^ — > C™'^" be the canonical projection which maps each 
point of 0,"^+^+^ on its first m + n coordinates. Then tt o Xj^ ,^ : W L,n x C" — >■ 
VL,n X C" is a finite bijective and birational morphism of affine varieties and 
therefore Al,„ : W^L,n x — Bin has the same property. This implies 
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'^{BL,n) = VL,n and that B^ n represents a rationally parameterized family of 
polynomial functions which extends the family represented by A^,™- 

Let Ki^n the function field over C of the absolutely irreducible variety V^.n 
and let Ri^n be the C- algebra of all rational functions 9 of Ki^n such that for 
any point u G Vi^n and any place ip : Ki^n — C U {00} whose valuation ring 
extends the local ring of the affine variety VL^n at the point u, the value (f{9) is 
finite and uniquely determined by u. Thus, for u G V^.n and : Ki^n — )■ C U 
{00} as above, we may associate to any polynomial f '■= J2 (^ii-i„x^i • • • ^n" ^ 
RlA^i, ■■ ■,Xn] the polynomial f{u, Xi, . . . x„) := Yl aii...i„{u)x\^ • • • := 

'^(«n...i,Ja^i' • ■_^n, which belongs to C[xi, . . . , 

Since : WL,n Vi^n is a finite, bijective and birational morphism 
of affine varieties, we conclude that RL,n contains the coordinate ring of the 
affine variety WL,n (see e.g. [H]). This implies that there exists a polynomial 

X\ , . . . Xn ] with the following property: for any 
u £ X G and y G C, the point {u,x,y) belongs to BL,n if and only 

if fL,n{u,x) = y holds. 

A branching-free output representation of the extended sample point 
query associated to the constructible set Ai^n can now easily be reahzed 
by any arithmetic circuit which first computes all monomial terms of the 
polynomial f^n and finally sums them up. 

Therefore, the given branching-parsimonious query evaluation procedure 
V produces on input consisting of the database D and the formula $l,„(S', P) 
a branching-free representation of the extended sample point query associ- 
ated to A^^n- This branching-free representation is realized by an essentially 
division- free single-output arithmetic circuit /3L^n with inputs Xi, . . . ,x„ and 
scalars 6[^''^\ . . . , belonging to Ri^n such that f3L,n computes at its out- 
put the polynomial fL,n e /2l,„[xi, . . . , x„]. 

Let t and £1 be new variables and let us now consider the polyno- 

mial gn '■= ^ ni<j<n(^* + -^i) defining the constructible object class 

:= {r Yl i^i + ^i) I Ai, . . . , A„ G C} 

l<j<n 

of n-variate complex polynomial functions. 

Observe that each element of r„ has nonscalar sequential time complexity 
at most n. 

Therefore the Zariski-closure r„ of the object class is contained in 
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Observe that is an absolutely irreducible, Q-definable affine variety. 
Since : Wn,n Vn,n is a finite morphism of irreducible affine varieties, 
we conclude that C„ := /in,n(r„) is an absolutely irreducible, Q-definable 
closed affine subvariety of Vn^n- 

Let En and L„ be the coordinate ring and the rational function field over 
C of the absolutely irreducible affine variety C„. 

Observe that we may identify En with C[(y^„(t, £i, 71), 5^„(t, 

il, ■ ■ ■ ,L,7m)] and Ln with C^gnitJi, . . . ,in,Jl); ■ ■ ■ i9n{t,il, . . . ,in,7m))- 

Therefore we may consider En as a C-subdomain of the polynomial ring 
C[t,£i, ... ,4] and L C-subfield of C(t, 4, . . . , 4)- 

The rational functions ^^"'"^ , ■ ■ ■ , Oi^'^^ of the affine variety Vn,n may 
be not defined on the subvariety C„. Nevertheless, since they belong to 
the C-algebra Rn,n, one verifies easily that there exist rational functions 

, • • • , o"l"^„ of the affine variety C„ satisfying the following condition: for 
any point w G C„ and any place i/j : L„ — )■ CU{oo} whose evaluation ring 
extends the local ring of Cn at the point u, the values of iIj at crj"^ , ■ ■ ■ , 
are given by ipic^i^^) = 0^'"'\u), . . . , ■0(cri")„) = 0'i^'^\u) and therefore finite 
and uniquely determined by u. 

In particular, the rational functions , • • • , G L„ are integral over 
the C-algebra En and hence contained in the polynomial ring C[t, ii, . . . , in] 
(see, e.g., [IS]). Therefore we may consider a["'\ . . . , as polynomials in 
the variables t and ii, . . . ,in (i.e. as elements of C[t, ii, . . . , in]). 

^From Qn = ^ ni<j<n(^« + -^j)' infer the identities 

Since the polynomials a["'\ . . . , ai^\ depend integrally from gn{t,ii, 
■ ■ ■ ,^n, Im), • • • , dnii, £i, . . . , 7^), wc dcducc uow casily that the polynomi- 
als {0,ii, . . . ,in), ■ ■ ■ , cri^]^ {0,£i, . . . ,in) do not depend on the variables 
ii, . . . ,in, i.e^ they belong to C. 

Let now be the division-free arithmetic circuit with scalars in C[t,£i, 
. . . ,in] obtained by replacing in the circuit f3n,n the scalars ^j"'"^ , . . . , 6*1"'"^ 
by the polynomials ci^i*,, . 

One verifies easily that the circuit /3„ computes the polynomial 

= i J]^ (^i + Xi) 

l<i<n 



LCi . . . . . . d^n 

«l,...,5„,ei,...,e„e{0,l} 
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Let vi, . . . , Vs„ „ be new variables. From the directed acyclic graph struc- 
ture of (3n (or I3n,n) ouc dcduccs immediately that for each {Si,...,5n) € 
{0,1}" there exists a polynomial G Q[vi, . . . ,Vs„„\ satisfying the 

condition • • • ' ^^i) = ■■■^n- Let Qn : C^'- ^ C^" the 

polynomial map defined by Qn := (^i, • • ■ , (^n) e {0, 1}"). 

Consider now an arbitrary integer 1 < p < 2" and let Ap^i : = 
p2'\...,Ap,„ := p2"-'^ ;= (Ap,i,...,Ap,„) and aj,"^ : C C^"-" and 
: C —7- C^" be the parameterized algebraic curves defined for r G C 

by 

:= (ai"Hr,A,),...,<)jr,A,)) 

and 

Pir\r) := (rA^;, . . . Aj- . . . , 5„) G {0, 1}") = (V; < j < 2"). 
Observe that the functional identity 

= g„oa(«) (3) 

is valid. 

Since the polynomials o"("^(0, £i, ...,£„),... , cri^)„(0, ^i, . . . , £„) do not de- 
pend on the variables £i, . . . , there exists a point a„ G C'^"-", independent 
on p, such that a^f^\o) = a„ holds. 

Let us denote the derivatives of ctp"'', and by and -DQn- 

Furthermore let Wp"^ := ^^^(0) G C''"-", let //„ : C"'" -> C^" be the C-linear 

map defined by := (DQ„)(a„) and observe that -^(0) = (/>'; < j < 2"-) 
holds. Applying the chain rule to ([3]), we infer the following identities: 

Since (p^)i<p<2",o<j<2" is a nonsingular Vandermonde matrix, we conclude 
now that the image of the C-linear map rjn '■ C'^" " — )■ C^" contains 2" linear 
independent points. Therefore rjn is surjective. This implies Sn,n > 2"-. 

Therefore the arithmetic circuit which represents the output pro- 
duced by the procedure V on input consisting of the database D and the 
formula $n,n; contains at least 2" scalars. 
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This implies that the nonscalar size of the circuit /3„^„ is at least 2 2 — n— 1. 
In conclusion, the procedure V requires 2^^"^ sequential time in order to 
produce the output /3„^„ on input consisting of the data base D and the size 



The main outcome of Theorem [T] and its proof can be paraphrased as 
follows: constraint database theory applied to quite natural computation tasks, 
as, e.g., branching-parsimonious interpolation of low complexity polynomials, 
leads necessarily to non-polynomial sequential time lower bounds. 

In view of the Pr, 7^ A^-Pr conjecture in the algorithmic model of Blum- 
Shub-Smale over the real and complex numbers, it seems unlikely that this 
worst case complexity behavior can be improved substantially if we drop 
some or all of our previously introduced requirements on queries and their 
output representations. Nevertheless we wish to stress that these require- 
ments constitute a fundamental technical ingredient for the argumentation 
in the proof of Theorem [H 

7 Conclusion and future research on the com- 
plexity of query evaluation 

In this paper, we have emphasized the importance of data structures and 
their effect on the complexity of quantifier elimination. 

However, the intrinsic inefficiency of quantifier-elimination procedures 
represents a bottle-neck for real-world implementations of constraint data- 
base systems. As we have argued, it is unlikely that constraint database sys- 
tems that are based on general purpose quantifier-elimination algorithms will 
ever become efficient. Also, restriction to work with linear data, as in most 
existing constraint database systems [IHl Part IV] , will also not lead to more 
efficiency. A promising direction is the study of a concept like the system 
degree, that has shown to be a fruitful notion for the complexity analysis of 
quantifier elimination in elementary geometry and has been implemented in 
the software package (polynomial equation solver "Kronecker", see [H]). In 
the context of query evaluation in constraint databases, the notion of system 
degree is still unsatisfactory since it is determined both by the query formula 
and the quantifier-free formulas describing the input database relations. It 
is a task for future constraint database research to develop a well-adapted 
complexity invariant in the spirit of the system degree in elimination theory. 




□ 
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Another direction of research is the study of query evaluation for first- 
order languages that capture certain genericity classes. For example, the 
first-order logic F0( between) has point variables rather than being based on 
real numbers and it captures the fragment of first-order logic over the reals 
that expresses queries that are invariant under affine transformations of the 
ambient space [12]. Although a more efficient complexity of query evalua- 
tion in this language cannot be expected, it is interesting to know whether 
languages such as F0( between) have quantifier elimination themselves (after 
an augmentation with suitable predicates). 
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Appendix A: Language and tools from alge- 
braic geometry. 

Let k be the field Q, R or C. We assume k to be "effective" with respect to 
arithmetic operations as addition/subtraction, multiphcation/division. Let 
k be an algebraically closed field containing k (in the sequel we shall call such 
a field an algebraic closure of k). 

We denote by N the set of natural numbers and by Z+ the set of non- 
negative integers. 

Fix n e Z"*" and let Xq, . . . , be indeterminates over k. We denote by 
A"^ := A.^{k) the n-dimensional affine space and by P" := P"(A;) the n- 
dimensional projective space over k. The spaces A" and P" are thought 
to be endowed with their respective Zariski topologies over k and with 
their respective sheaves of fc-rational functions with values in k. Thus the 
points of A"^ are elements (xi, . . . ,a;„) of k and the points of P"^ are (non 
uniquely) represented by non-zero elements (xq, . . . , of /c and denoted 
by (xq : ••• : The indeterminates Xi, . . . ,Xn are considered as the 

coordinate functions of the affine space A". The coordinate ring (of polyno- 
mial functions) of A" is identified with the polynomial ring k[Xi, . . . 
Similarly we consider the (graded) polynomial ring k[Xo, . . . , X„] as the pro- 
jective coordinate ring of P". Consequently we represent rational func- 
tions of P" as quotients of homogeneous polynomials of equal degree be- 
longing to k[XQ, . . . , Xn]- Let Fi, . . . ,Fs be polynomials which belong to 
k[Xi, . . . , Xn] or are homogeneous and belong to k[Xo, . . . , X„]. We denote 
by {Fi = 0, . . . , Fs = 0} or V{Fi, . . . , Fg) the algebraic set of common ze- 
roes of the polynomials Fi, . . . , Fs in A" and P" respectively. We consider 
the set V := {Fi = 0, . . . , = 0} as (Zariski-)closed (affine or projective) 
subvariety of its ambient space A" or P" and call V the affine or projective 
variety defined by the polynomials Fi, . . . , F^. We think the variety V to be 
equipped with the induced Zariski topology and its sheaf of rational func- 
tions. The irreducible components of V are defined with respect to its Zariski 
topology over k . We call V irreducible if V contains a single irreducible com- 
ponent and equidimensional if all its irreducible components have the same 
dimension. The dimension dim V of the variety V is defined as the maximal 
dimension of all its irreducible components. If V is equidimensional we de- 
fine its (geometric) degree as the number of points arising when we intersect 
V with dim V many generic (affine) linear hyperplanes of its ambient space 



46 



A" or P". For an arbitrary closed variety V with irreducible components 
Ci, . . . ,Ct we define its degree as deg V := degCi + ■ ■ ■ + degCf. With this 
definition of degree the intersection of two closed subvarieties V and W of 
the same ambient space satisfies the Bezout inequality 

degVnW < deg V deg W 

We denote by k[V] the affine or (graded) projective coordinate ring of 
the variety V. If V is irreducible we denote by k{V) its field of rational 
functions. In case that V is a. closed subvariety of the affine space A" we 
consider the elements of k[V] as A;- valued functions mapping V into k. The 
restrictions of the projections Xi, . . . ,X„ to V generate the coordinate ring 
k[V] over k and are called the coordinate functions of V. The data of n 
coordinate functions of V fixes an embedding of V into the affine space A". 
Morphisms between affine and projective varieties are induced by polynomial 
maps between their ambient spaces which are supposed to be homogeneous 
if the source and target variety is projective. 

Replacing the ground field k by its algebraic closure k, we may apply all 
this terminology again. In this sense we shall speak about the Zariski topolo- 
gies and coordinate rings over k and sheaves of /s-rational functions. In this 
more general context varieties are defined by polynomials with coefficients in 
k. If we want to stress that a particular variety V is defined by polynomials 
with coefficients in the ground field k, we shall say that V is k-definable or 
k-constructible. The same terminology is applied to any set determined by a 
(finite) boolean combination of fc-definable closed subvarieties of A" or P". 
By a constructible set we mean simply a /c-constructible one. Constructible 
and fc-constructible sets are always thought to be equipped with their corre- 
sponding Zariski topology. In case of := Q and := C we shall sometimes 
also consider the euclidean (i.e. strong") topology of A" and P" and their 
constructible subsets. 

The rest of our terminology and notation of algebraic geometry and com- 
mutative algebra is standard and can be found in [IS], [22], and in [T]. 
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